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Abstract 

We report theoretical and experimental studies of the longitudinal electron spin and orbital 

£ ■ 

relaxation time of interstitial Li donors in 28 Si. We predict that despite the near-degeneracy of 

the ground-state manifold the spin relaxation times are extremely long for the temperatures below 

i i ' 0.3 K. Lhis prediction is based on a new finding of the chiral symmetry of the donor states, which 

t— I ■ presists in the presence of random strains and magnetic fields parallel to one of the cubic axes. 

> ; 

Experimentally observed kinetics of magnetization reversal at 2.1 K and 4.5 K are in a very close 
agreement with the theory. Lo explain these kinetics we introduced a new mechanism of spin 
decoherence based on a combination of a small off-site displacement of the Li atom and an umklapp 
phonon process. Both these factors weakly break chiral symmetry and enable the long-term spin 



relaxation. 



PACS numbers: 03.67.Lx, 76.30.-v, 71.55. Cn, 71.55.Cn 



I. INTORDUCTION 

Recent progress in building devices for quantum information processing (QIP) has lead to 
renewed interest in studying the relaxation times of electron spins bound to impurities and 
in quantum dots to be used for encoding quantum information^ 2 -^^. Electron spins bound 
to shallow donors in Si are characterized by particularly long relaxation times, since the 
spin orbit interaction is weak and isotopically enriched 28 Si is available^ 7 -. However a weak 
interdonor spin coupling represents a challenge for the controlled interaction of quantum 
bits. 

The interstitial lithium donor in silicon is unique amongst shallow donors because it has 
an inverted level structure and near degeneracy of the Is ground state manifold^. The 
latter gives rise to a strong long-range elastic-dipole coupling of the neutral Li donor orbital 
states^, similar to that obtained earlier for acceptors in Si^ but absent in the case of other 
donors. In^ a quantum computing architecture was proposed based on the stress-defined 
orbital qubits and the above interaction was used to enable 2-qubit gates. At the same 
time the near degeneracy of Is manifold in Si:Li also gives rise to a nontrivial coupling 
between the spin and orbital degrees of freedom first investigated by Watkins and Ham^. 
One can expect that the inverted level structure and near degeneracy in Si:Li will lead to 
a very short electron spin relaxation times which can also strongly affect orbital relaxation 
and be detrimental for QIP. Somewhat surprisingly, in view of extensive studies of spin 
relaxation of other shallow donors in silicon since the early sixties, both experimental^ 7 -^^ 
and theoretical 1 ^, the question about Li spin relaxation has not been addressed. 

In this work we propose a new relaxation mechanism for coupled donor electron spin and 
orbital degrees of freedom. It is based on (i) our new finding of the chiral symmetry of 
the donor states that persists in prsence of the random starin and magnetic field of certain 
orientations and (ii) an interplay between the umklapp phonon processes and a short-range 
electric dipole interaction that only weakly break this symmetry and enable spin relaxation. 
As a result Li electron spin and orbital relaxation times become very long in the presence 
of modest external stress (~ 10 7 Pa) at temperatures below 0.3 K. We also report pulsed 
ESR measurements of the longitudinal spin relaxation time Ti in 28 Si:Li at temperatures 
1.8 to 4.5 K. We find that in the absence of external stress these times are shorter than our 
instrumental resolution (< 100 ns) at temperatures as low as 1.8 K. This rapid relaxation is 



caused by the degeneracy in the ground state manifold. However, a modest applied stress lifts 
the degeneracy, and leads to a much longer relaxation, as will be discussed below. The results 
of the theoretical analysis of the selection rules and spin-orbital relaxation mechanisms are 
in excellent quantitative agrement with the experiment as shown in Fig[U 

The ground state Is manifold of an interstitial lithium donor in silicon can be described 
by the Hamiltonian in a twelve dimensional Hilbert space formed by a tensor product of 
electron spin states |x±) and six Is orbitals corresponding to the six conduction band minima 

H = H<® + H d , H® = H + H B + H st + H so . (1) 

The term H Q is the spin-diagonal, effective mass orbital Is Hamiltonian also describing the 
valley-orbit splitting due to the short range tetrahedral potential. The rest of individual 
terms in H to be considered below describe the coupling of the donor electron to magnetic 
field (H B ), long- wavelength strain field (H st ) , spin-orbit interaction (H so ), and a small 
off-center displacement of Li core from the symmetry site (Hd). 

For each of the two spin orientations the six eigenstates of H are linear superpositions 
of the six silicon valley orbitals, ^(r) = J2j=±i a jFj(. r ) u j( r )- Here Fj(r) and Uj(r) are, 
respectively, the hydrogenic envelope and the Bloch function for the j-th conduction band 
minima k, = k Xj (xj=±x,±y or ±z). Eigenstate ^(r) belongs to the irreducible rep- 
resentation /i of the tetrahedral group characterized by the valley-orbit coefficients a^. In 
contrast to the Group V substitutional donors the ground state of H is five-fold degenerate 
for each spin orientation. It is composed of a triplet of states T 2x ,T2 y ,T 22 which are anti- 
symmetric combinations of just two opposite valleys, a^ = —octj, and a doublet Eg, E e that 
are symmetric combinations of valleys, af* = a^_ •, with a second-order zero at the origin. A 
fully symmetric singlet state A 1 lies at A = 1.76 meV above the ls(E + T 2 ) orbital ground 
stated, so that matrix elements (V^j-HolVy) = ^o^^Ai^^'Ai^ where I is 2x2 unit matrix 
operating in spin subspace. 

II. PARITY-SELECTION RULES FOR ACOUSTIC-PHONON TRANSITIONS 

Electron-phonon interaction H^ within ls(E+T2) manifold can be described in terms 
of deformation potentials because the states there decay via emission of long-wavelength 
acoustic phonons with g~10 5 cm _1 . A matrix element between a pair of donor states ip^ and 





—* — '■ ■* ,_ ~ 


1.5 




A 


^^ 




(a) - 


> 1 








— 


S, 




\ 


& 0.5 

■— 




* 2x>hy'~e^e 


a 






w 


"^=^ E e ~ 


-0.5 


- . 2z . r^ 







10 







(b) A,/ 


3 






V 


2/ 




il_. 


' > 


' 



T T 



"6i 



T 



2zt 



r 



o 



F 100 (10 Pa) 



2zi 



8-10 6 

■(-' 

•a 

t 4-10 6 

I ° 


|'|| 


#^ T = 4.5 K 

/ , i , i 


] 


Temperature (K' 
L0 0.25 O.l 


- 


r2 


io 6 

^ 2 

£10 

H -2 

10 


i i i 


- 


O 

S "4-10 

N 

•iH 
-M 

fl 

«3 -8-10 


- 2y^ l - 

/A 2 =0.llmeV 

~ 3/ 

_/A 3 =0.54meV 

i/ i i 


I 


l/T — *- 

i 





10 15 

Time (|isec) 



FIG. 1: (a) Theoretical stress dependence of Si:Li Is energy levels at B = 0. (b) Energy levels at 
finite stress and B / 0. Solid lines 2,3 indicate the relatively fast one-phonon transitions connecting 
the states with the same chirality. Dashed line 1 indicates the slow spin-flip transition between the 
"left" and "right" states; (c) Comparison of experimental and theoretical kinetics of magnetization 
recovery at 2.1 K and 4.5 K. The theoretical curves (thin solid lines) correspond to Vg = V sin.9 = 7x 
10~ 4 meV and Fioo = 9 x 10 6 Pa. The lower inset displays the theoretically calculated temperature 
dependence of the spin lifetime r = 1/A m i n , an inverse of the lowest eigenvalue of relaxation operator 
£, Eq. (fTOl) , It is characterized by two activation energies A2 and A3 (the corresponding activation 
processes are shown in Fig. CUb) with dashed lines) 
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■0' of Hq U with phonon wavevector q and polarization v obeys parity selection rules^ and 
can be written in the form 

3 

<W|tf q ,|Vy > = / £ M ^ A L,' + W**jw?\> ( 2 ) 

where A^ , = otfot^j ± a^a^-. Here M^„ is the energy shift of valleys k±j due to the 
phonon mode q with polarization ia In can be expressed in terms of the deformation 
potential constants and angular functions^. A + and A~ are zero for the states /i, //' with 
the same and opposite parities, respectively. In the latter case the transition amplitude ([11) 
contains a very small quantity qr u ~ 10~ 4 where r u = J F^(r)zsm(K u z)dr is determined by 
the overlap in k-space of the envelope functions localized near the different conduction band 
minima separated by umklapp vector k u ~ 0.3 x 10 8 cm" 1 — . 

Matrix elements {ip^\H st \ip^i) = 2„J V =1 At «,., depend on a shear deformation poten- 
tial E u that couples the long-wavelength strain tensor components Ujj to the donor states^. 
For most experimental situations the strain is caused by a superposition of an external uni- 
axial and much weaker random biaxial stress that completely removes the orbital degeneracy. 
A compressive uniaxial (001) stress F oi only partially splits the ls(E + T 2 ) manifold*^ as 
shown in Fig. [TJ The ground state is T 2z and the first excited state is a superposition of 
Eg and A x orbitals with dominant Eg character. The third level is triple-degenerate, corre- 
sponding to states T 2x , T 2y , E e as shown in Fig. [Q For F 00 i = 9xl0 6 Pa the splitting of T 2z 
and Eg levels is ~ 0.2 meV; much larger than the random-stress induced level broadening 
in 28 Si — <1&. We note that due to the symmetry of an interaction involving deformation 
potentials^ an arbitrary biaxial stress mixes only even-parity states, Ai,Eg and E e and 
therefore preserves the above parity-selection rules. 

Matrix elements of the Zeeman Hamiltonian equal 

3 

(Vv|F B |Vy> = Y.S^bB.S^^, + eA+^,), (3) 

where \ib is the Bohr magneton, Sj and Bj are, respectively, the 3 cartesian components 
of the electron spin operator and magnetic field. Also e = (g\\ — g±)/g±, and g\\,g± are the 
components of the g-tensor— . 

There exists a small spin-orbit interaction in the ground state of Li donor induced by 
the host silicon lattice^. It is typically smaller then the strain and Zeeman terms. The 



spin-orbit Hamiltonian H so reads^: 

H so = -iS z (X'\T 2z )(E £ \ + X\T 2x )(T 2y \) 
-iS x (X'\E + )(T 2x \ + X\T 2y }(T 2z \) 
-iS y (X'\E.){T 2y \ + X\T 2z ){T 2x \) + h.c, (4) 

where \E±) = \ (\E E ) ± V3\E e )), and small constants A ~ 0.003 meV and A' = 0.007 meV 
determine the strength of the spin-orbit coupling within the triplet and between the triplet 
and singlet, respectively^. 

III. CHIRAL SYMMETRY 

Spin-orbit Hamiltonian H so mixes the orbital states of opposite parities and, as one 
might expect, it violates the above parity selection rules for the transition rates. However, 
we show below that the time-reversal symmetry in the absence of a magnetic field leads 
to a factorization of the total Hamiltonian and emergence of the new selection rules that 
persist for finite magnetic fields of certain orientations. We assume that the magnetic field 
B = Bz. Then the eigenstates of the Hamiltonian H^ form two non-overlapping subspaces 
with an element in each subspace being a linear superposition of the following 6 basis vectors 

\Xo)\^»)- 

l^ 0) ) = \Xa) (c z \T 2z ) + c e \E e ) + c e \E £ ) + c Al \Ai)) 

+ \x-a) (c x \T 2x ) + Cy\T 2y )) , a = ±1, (5) 

where & z \Xa) — ^\Xa)> °~z is the Pauli matrix, and a is a subspace index. The eignstates 
|$+ ) and their time-reversed counterparts |$_ ) are also the eigenstates of the "chirality" 
operator Z = R z (ir)a z , where R z {n) is the operator of rotation through the angle n around 
2-axis. Indeed, the operator Z commutes with the Hamiltonian H^ and Z\$>±) = ±|$± ). 
In the absence of the magnetic field pairs of similar eigenstates with opposite chiralities 
form Kramers doublets with the same energy. A finite magnetic field B lifts the Kramers 
degeneracy. However for Hb with B = Bz the chirality still remains a good quantum number 
and the subspaces do not mix despite an arbitrary biaxial strain present in H^\ Similar 
subspaces can be constructed for B = Bx and B = By. 



Watkins and Ham postulated a small displacement d of the Li atom from a tetrahedral 
interstitial site^, preferentially in the direction of the uniaxial stress F || d to generate the 
correct number of lines in the Si:Li ESR spectrum. It gives rise to a short-range electric 
dipole potential H d (PQ) with the strength V and matrix elements 
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<Wl#d|Vy > = 2iVI ^ Aj^dj/d, dj = d- Xj, (6) 

connecting states of opposite parities and the same spin orientation. Using ((5|), (J6j) and 
coefficients atj one gets 

($(o)|#|$W) = for B || d || Xj, (7) 

where Xj = z, y or x. However, we will assume a small random misalignment, \xj x d/d\ = 
| sin 9 1 <C 1 (perhaps, due to a random biaxial stress in addition to a strong uniaxial stress 
F=Fx j ). It will lead to a slight mixing of the "left" and "right" states ($S 0) |# |$L°i) ex Vg = 
Usinfl. 

IV. TRANSITION RATES 

Consider B||F||xj (xj = z,y or x). Then the matrix elements of H qL/ between a pair of 
exact eigenstates, H§ n = £ n $ n , that are predominantly "left" and "right" have the form: 
(<& n \Hqy\<& m ) — C^ u qr u +C^ u £ J* £ . The first term is due to an umklapp process ( gr u ~10~ 4 ) 
and the second term is due to a small misalignment Vq=V sva.9 of d and Xj. Both terms are 
small and acoustic phonon transitions between the states of opposite chiralities are suppressed. 
As shown below, under the appropriate conditions each of the chiral subspaces can be made 
nearly "decoherence free". 

An important application of the above rule is a suppression of the spin-flip transitions 
between the states with dominant characters |x±i) <S> |Vv) that differ only in spin orientation. 
Those states always have opposite chiralities We shall label each eigenstates after its domi- 
nant character and consider B = Bz, F = F 00 iz. If the compressive stress F 00 i is sufficiently 
large the Zeeman splitting g^HsB becomes smaller than a stress-induced splitting between 
Eg and E T2z (see Fig. [2]). Then at low temperature the state |T 2z |) decays directly into the 
state |T 22 |) of the opposite chirality. There is a competition between umklapp and dipole- 
induced processes. For moderate uniaxial stress and magnetic field the dipole contribution 
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FIG. 2: (Color online.) (a) Calculation of the first four energy levels of Si:Li donor in a magnetic 
field B = 1T; (b) calculation of the lifetimes of the first two excited states versus compressive 001 
stress at T=0.1 K, B = 1 T, and Vg = V s'rnO = 0. Solid lines — B || z and dashed lines — B || x. 
Inset shows magnified region of the avoided crossing for the magnetic field B JL z; (c) calculation 
of the lifetimes of the first two excited states versus compressive 001 stress at T=0.1 K, B = 1 T, 
and Vf> = V sin9 = 4.8 x 10~ 4 meV. Solid lines — B II z and dashed lines — B II x. 



to the lifetime tt 27 t dominates 



1 (\\'V d sm6 d y 4lv 3 ( (x 



<^rt\'^P^\ ^ r- + V » . (8) 



T *X~ M \ A 3 J x 8 \\3 
Here v = g±fiBB/A and x = E u (Su — Sii) -Fooi/6 A where Sn, Su are elastic compliance 
moduli. Prefactor in (JSj) corresponds to the lifetime Ta x ~1 psec of the states A\ that decay 
within the same subspaces to Eg and E e . 

For greater stress and/or magnetic field the umklapp process dominates the decay, giving 

/AA'AoaiiV 7v 5 



r T 2zT - T A 1 iu ( " 



hu± J 3(x 2 — v 2 ) 2 



Here a\\ is the longitudinal Bohr radius, u± is the transverse sound velocity, and £ u ~ 5- 1CT 5 
depends on the overlap factor r u given above via {r u /2a\\) 2 = £ u — £,/ . 

We calculated transition rates W n ^ m between each pair of the 12 eigenstates |$ n ) by 
numerical diagonalization of H (pp) and employing the electron-phonon interaction ([2]). The 
stress dependence of the lifetimes of the first two excited states of Si:Li donors at low 
temperatures is shown in Fig. [2]for B || z and B || x. These states can be classified as either 
orbital (predominantly Eg^) or spin (predominantly T 2z ^) excitations. If the stress-induced 
splitting (SIS) is smaller than the Zeeman splitting (ZS) the first and second excited states 
correspond to the orbital and spin excitations respectively. As soon as SIS overcomes ZS 
the two states switch their characters. 

The orbital and spin excitations have distinctly different stress dependencies of their 
lifetimes. The lifetime of the orbital exitation is decreasing with the stress increase while 
the spin excitation behaves quite oppositely. The stress increase leads to a large separation 
between E$i and T 2z i states. Consequently the probability of E$i — > T 2z i phonon decay 
grows as the phonon density of states or faster. On the contrary, the stress does not affect 
the Zeeman doublet T 2z -\ — T 2z ^ but rather separates it from other states of the manifold 
making the spin excitation T 2z ^ more stable. The lifetime decrease of the orbital excitation 
is more pronounced for an ideal position of the Li-donor, i.e. for Vg = (see Fig. E](b)) and 
for B || x. In this case the orbital excitation and the ground state have different chiralities 
and the former may decay through the umklapp process only. This adds additional stress 
dependence due to the matrix element of the umklapp process. The magnitude and the stress 
dependence of -E^-lifetime shown in Fig. [2](b) (dashed line) is very similar to our previous 
result^. The dipole term Vg reduces both the magnitude and the stress dependence of 
-Eej-lifetime, which is weakly decreasing with stress near the value of 10" 4 s (see Fig. [2(c)). 
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V. EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

We conducted a pulsed ESR study of the Li donor spin relaxation under stress. An 
isotopically-enriched 28 Si crystal (800 ppm residual 29 Si) was 7 Li doped by implanting with 
a variable energy, up to 360 keV, and a total dose 6T0 12 of 7 Li/cm 2 followed by a 30 min 
anneal at 800 °C to ensure a homogeneous distribution of the lithium. The estimated Li 
concentration is ~ 10 14 cm -3 . Uniform, (110), in-palne tensile stress was applied to the Si by 
gluing fused silica slides (1 mm thick) to the two (001) faces with Apiezon N vacuum grease. 
This "sandwich" was then frozen to liquid helium temperatures. The thermal expansion 
coefficient of fused silica is substantially smaller than that of silicon in the temperature 
range 4-300 K— , and therefore freezing the sandwich results in an uniform in-plane tensile 
strain, equivalent to the action of a compressive uniaxial stress F oi along the (001) axis of 
the silicon crystal (perpendicular to the faces). The spin resonance experiments were done 
at X-band (~0.34 T). Under F oi the ESR spectrum is a single line as described by Watkins 
and Harn^ though our linewidth (8 //T) is considerably narrower than in their natural Si, 
since inhomogeneous broadening from the 29 Si is eliminated in the isotopically enriched 
28 Si^. The electron magnetization decay measured at 2.1 and 4.5 K with a conventional 
inversion recovery pulse sequence^ 1 - is shown in Fig. [TJ 

The stress F om in the sandwiched samples was estimated by performing Fourier transform 
infrared spectroscopy (FTIR) with a Bruker IFS 66v/S spectrometer (resolution, 0.14 cnU 1 ) 
on a similarly prepared stressed samples of Li-doped natural float zone silicon (a hole with 
a diameter 3 mm in the silica slides allowed FIR transmission through the stressed sample). 
The stress-induced splitting 0.202 meV was observed between the ls(E + T 2 ) — > 2p (—) line 
at 21.287 meV and ls(E + T 2 ) — > 2p line at 21.489 meV. Using this splitting we estimate 
F oi — 9.0 x 10 6 Pa based on the previous measurements with calibrated stresses 9 - and 
expressions for the stress-induced line shifts^ 9 -. 

To explain the magnetization recovery data in Fig.[U(b) we studied the multilevel kinetics, 

12 

Pn / *~-nmPmi [^ ) 

m=l 

where 

12 



fc=l 
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is a relaxation operator and p n = ($ n |p|$ n ) are individual state populations. We used the 
stress value Fioo = 9x 10 6 Pa estimated from the above FTIR experiment. The only unknown 
parameter in the expressions for the decay rates W nm is Vq, or rather V e 2 . We were able to fit 
the two experimental time traces at 2.1 K and 4.5 K using Vg = 7x 10~ 4 meV (~ 0.25A) as 
shown in Fig. [U This value corresponds to a statistical average {{V e 2 )) 1 ^ 2 over all Li donors. 
We note the multi-exponential character of the decay at 4.5 K. The temperature dependence 
of the long-time magnetization kinetics is determined by that of the smallest eigenvalue of 
the relaxation operator A m i n (T). The magnetization decay proceeds via thermal activation 
to rapidly decaying excited states as shown in Fig. [Tj(b). For T < 0.1 K the spin life-time r 
saturates at large value consistent with those shown in Fig. EJ 

VI. CONCLUSIONS 

To conclude, we introduced a new spin-orbital relaxation mechanism for Li donors in Si 
and predicted theoretically that the longitudinal spin relaxation 7\ time of a Li donor in Si 
under the uniaxial compressive stress and magnetic field aligned along the x,y or z axis can 
be very large (hours) for T < 0.3 K despite the inverted level structure of interstitial Li donor 
as compared to Group V substitutional donors in silicon. The orbital relaxation time is 0.1 
msec at that temperature. By virtue of the chirality selection rules for the electron-phonon 
interaction, the donor electron spin can decay only via inter-valley umklapp processes or a 
very weak direct phonon process caused by a small Li off-site displacement misalignment. 
At higher temperatures T\ is reduced to a /zs range due to thermal population of the lower 
excited states which in turn promotes fast transitions within the same chiral subspace allowed 
by the parity-selection rules. We obtained a very close (within few per cent) agreement with 
pulsed ESR measurements of T\ times at T = 2.1 K (12 ps) and at T = 4.5 K (1.5 /zs) that 
corroborates our prediction of extremely long-lived spin states of Li donors in Si at lower 
temperatures. Since the chiralty is a universal property of the silicon host ur study can be 
applied to substitutional donors in Si with similar near-degeneracy of the ls(E + T 2 ) coupled 
to the ground state via stray electric fields and umklapp processes. Our results of the very 
long Li spin lifetime can also stimulate the analysis of the possibility of Li spin-based QIP 
applications with direct long-range elastic-dipole interaction between the Li donor spins 
similar to that considered ir*^ 1 ^. 
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